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Abstract: A high significance observation of two muon signals, the shadow of the sun and moon, have 
been seen by the 5.4 kt MINOS Far Detector, at a depth of 2070 mwe. The distribution of angular separa- 
tion of muons near the moon was well described by a Gaussian, which was used to determine the angular 
resolution (0.34° ± 0.07°) and pointing (0.3° ± 0.05°) of the detector. 
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Introduction & Motivation 

The MINOS Far Detector is a magnetized scintil- 
lator and steel calorimeter, located in the Soudan 
Mine in northern Minnesota, USA, at a depth of 
720 m. While primary function of the Far Detec- 
tor is to detect neutrinos from Fermilab's z/^ beam, 
the great depth and wide acceptance of the detector 
combined with the flat overburden of the Soudan 
site allow it to serve as a cosmic -ray muon detec- 
tor as well. The detector is composed of 486 8 m 
octagonal planes 2.54 cm thick, spaced 5.96 cm. 
This 5.4 kt detector is 30 m long and has a total 
aperture of 6.94 x 10'^ cm^ sr. MINOS is the first 
underground experiment able to discriminate pos- 
itively charged particles from those that are neg- 
atively charged for the purpose of CPT violation 
investigations, but this feature also allows inde- 
pendent study of positively and negatively charged 
cosmic rays. 

Optical telescopes use a standard catalogue of stars 
to establish the resolution and pointing reliabil- 
ity of a new instrument. This is not possible 
for a cosmic ray detector, as there are no cos- 
mic ray sources available for calibration. There 
are two well observed phenomena in the otherwise 
isotropic cosmic ray sky, though they are deficits, 
not sources. The sun and the moon provide a 
means to study the resolution and pointing of a cos- 
mic ray detector because they absorb incident cos- 
mic rays, causing deficits from their respective lo- 



cation. These signals allow a measure of phenom- 
ena associated with cosmic ray propagation and in- 
teraction resulting from geomagnetic fields, inter- 
planetary magnetic fields, multiple Coulomb scat- 
tering, etc. These extra-terrestrial objects have the 
same 0.5° diameter as viewed from Earth, though 
the sun shadow is more difficult to observe because 
is much farther away and has its own magnetic 
field that deflects the charged cosmic rays. The 
cosmic ray shadow of the moon has been measured 
by air shower arrays (CYGNUS [1], CAS A [2] 
, Tibet [3] ), as well as underground detectors 
(Soudan 2 [4], MACRO [5, 6], L3h-C [7]). 

Data 

This analysis encompassed events recorded over 
1339 days, from 1 August 2003 - 31 March 2007, 
for a total of 1194 live-days, and includes 51.41 
million cosmic ray induced muon tracks. Cos- 
mic ray muons were triggered by recording hits 
on 4/5 planes or exceeding a pulse-height thresh- 
old and were written to a temporary disk at Soudan 
and later sent to Fermilab for reconstruction. Sev- 
eral cuts were required to ensure only well recon- 
structed tracks were included in the analysis. Pre- 
analysis and run cuts include: failure of demulti- 
plexing figure of merit and data taking quality [8]. 
Analysis cuts include: track length less than 2 m, 
number of planes less than 20, Xreco/''^4f > 1-0, 
and either track vertex or end point outside of the 
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fiducial volume of the detector. A total of 30.52 
million events survived these cuts for the combined 
sample. 

The background for this analysis was calculated 
using a simple Monte Carlo simulation that ex- 
ploits two key features of the muons induced by 
cosmic ray primaries: the time between consecu- 
tive cosmic ray arrivals follows a well known dis- 
tribution [9] and the cosmic ray sky is isotropic, 
Thus, a bootstrap method that independently 
chooses the arrival time and location in space ef- 
ficiently simulates a cosmic ray muon. This sim- 
ulation chose a muon out of the known distribu- 
tion of events in the detector (in horizon coordi- 
nates), paired it with a random time chosen from 
the known time distribution, and found the muon's 
location in celestial coordinates. This was done 
for every muon to create one background sample; 
500 background samples were simulated for a high 
statistics background distribution. 

Combined Muon Calibration 

A further cut, excluding tracks with p < 
30 GeV/c, was required to reject muons that were 
greatly affected by energy dependent processes. 
The final data set included 20.17 milhon muons. 
To find the one dimensional space angle separation 
from the moon and the sun for each muon track, a 
hst of the celestial body's location in celestial coor- 
dinates in one hour increments was obtained from 
the JPL HORIZONS [10] ephemeris database for 
May 1, 2003 until May 1, 2007. The ephemeris 
data was referenced to the location of the detec- 
tor in latitude, longitude, and distance below sea 
level. A function was written that interpolated the 
celestial body's location at the time a particular 
muon is seen in the detector Then the particu- 
lar muon's arrival direction was compared to the 
celestial body's location at that time. The one di- 
mensional space angle separation was found using 
the Haversine formula, published in Sky and Tele- 
scope [?]. The muons were binned in 56i„ =0.10° 
increments, and since radial distance from the cen- 
ter of the celestial body is measured over a two 
dimensional projection, the bin solid angle of bin 
(i) increases when moving out from the center as 
Afli = {2i - 1) 

* SbinT^- Weighting the number 



of events in each bin by the reciprocal of the area 
resulted in the distribution Ni/Afl^, the differen- 
tial muon density. The distribution is shown in 
Fig. 1 (1); from the sun, Fig. 1 (r), with statistical 
error bars displayed. There is a very clear deviation 
from a flat distribution in both plots as — > 0, and 
that deviation is attributed to muons blocked by the 
moon and sun, respectively. The background was 
calculated using the method described in Sec. 2. 
As expected, the backgrounds are nearly flat, with 
a fit value of x^/nd/ = 0.12/39 for the moon and 
X^/ndf = 0.27/39 for the sun . This Monte Carlo 
is consistent with the premise of a sourceless cos- 
mic ray sky. The significance of the shadow can be 
found by fitting to the distribution a function of the 
form [4]: 

^=A(l-(i?ya^)e-^V2.^) (1) 

where A is the average differential muon flux, a ac- 
counts for smearing from detector resolution, mul- 
tiple coulomb scattering and geomagnetic deflec- 
tion, and Rm = 0.26°, the angular radius of the 
moon. A fit to eq. 1 yields /ndf ~ 37.9/38, an 
improvement of 16.4 over the liner fit (x^/ndf = 
54.3/39), with parameters A = 483.9 ± 3.1 and 
cr = 0.34° ± 0.07°. The change in over 38 
degrees of freedom for the moon corresponds to 
a 10~4 chance probabihty. The improvement of 
X^ probability for the sun corresponds to a 10~^ 
chance probability. These results are summarized 
in Table 1. 

Since the tracks of dimuon events, muons that are 
induced by the same cosmic ray primary, are nearly 
parallel when they are created, they can be used to 
find the MINOS Point Spread Function. To quan- 
tify the absolute pointing of the Far Detector, the 
MINOS Point Spread Function will be used to find 
the two dimensional contours of the most signifi- 
cant muon deficit caused by the moon. This anal- 
ysis is not complete at the time of this writing, 
however, so in its stead a simple approximation 
was used on the one dimensional moon shadow 
to place an quantify the pointing of the detector. 
A significant shadow pointing to the apparent lo- 
cation of the moon was found using a one di- 
mensional space angle separation, a convolution of 
ARA, ADec. Since the moon's radius is 0.26°, 
we can approximate the absolute pointing of the 
detector as 0.3° ± 0.05°, with the error given by 
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Figure 1: The differential muon flux with respect to displacement from the moon's (1), and sun's (r) location, 
binned in 0.1°. The dashed curve is the calculated background, while the solid curve is the best fit from eq. 1. 



one half the bin width in theta. This analysis as- 
sumes that the pointing is reliable to begin with, so 
more precise analysis is still to be performed. 

Charge-Separated Analysis 

For the charge separated sample, a further cut 
was required to exclude events with low confi- 
dence charge sign determination. The curvature of 
the track is used to determine the momentum and 
charge of the particle, so this cut was charge over 
momentum divided by the error in the determina- 
tion of charge over momentum (-2/2- > 2.2) [?]. 

Cg/p 

There are 2.7 million positively charged and 1.9 
million negatively charged muon in this sample. 
An analysis similar to what is described in Sec. 3 
was performed on the charge separated muon sam- 
ple. The charge separated moon shadow plots can 
be seen in Fig. 2 (1), and the sun shadow plots 
are shown in Fig. 2 (r). The Gaussian hypothesis 
gives very little improvement over the flat hypoth- 
esis, suggesting that the shadow of the moon or sun 
has not been seen for the charge separated sample. 
The strong cut reduces the statistics by over 
four times, from 20 to 4.6 million, eliminating an 
observable deficit in the direction of either sun or 
moon. This cut does give a charge ratio of 1.3 for 
both sun and moon distribution, which is consis- 
tent with the published MINOS result. This anal- 
ysis has provided a flat distribution in the region 
approaching the moon, so given more statistics a 
significant shadow should be observed. 



Conclusions 

Using 20.17 million muons accumulated over 1 194 
live-days, the MINOS Far Detector has observed 
the cosmic ray shadow of the moon with a high sig- 
nificance. Despite the inherent fragility of the one- 
dimensional moon shadow measurement (there 
were only nine events in the bin nearest to the 
moon), the null hypothesis of the muon deficit in 
the area near to the apparent location of the moon 
has a probability of 10^4. The cosmic ray shadow 
of the sun over the same time period has a chance 
probability of lO^'^The shadow of the moon was 
used to approximate both the effective angular res- 
olution of the detector, 0.34° ± 0.07°, and the ab- 
solute pointing of the detector, 0.30° ± 0.05°. In 
agreement with the expectation, no significant dif- 
ference was found in the shadowing effects for ei- 
ther population, save for the charge ratio. 
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Figure 2: The differential muon flux with respect to displacement from the moon's location (1), and sun's location 
(r), binned in 0.1° increments for fi'^ (open circles) and (open triangles). The solid curve is the best fit 
from eq. 1 for ; the dashed curve is for ii~ . The fit results are in Table 1. 
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